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Abstract 

In a recent experiment, Ma et al. measured the temperature dependence 
of the gap anisotropy of oxygen-annealed Bi2Sr2CaCu208 +a; . Their measure- 
ments were taken along the two directions T — M and T — X. They found that 
the gap along both directions is nonzero at low temperatures and that the 
ratio is strongly temperature dependent. We show, using Ginzburg-Landau 
theory, that this behavior can be obtained if one assumes the existence of 
s-wave and d-wave components for the order parameter. Our theory predicts 
orthorhombic anisotropy in the gap and anomalous behavior for the electronic 
specific heat below T c . 

PACS numbers:74.20.Mn,74.25.Nf,74.72.Bk,74.25.Bt 
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The most central issue in high-T c superconductivity at the present time is the nature 
of the order parameter. Many experiments have attempted to clarify the form and symme- 
try of the order parameter, with the most conclusive being of the phase-sensitive tunneling 
variety These indicate that the gap function changes sign when rotating by tt/2 in 

momentum space. This suggests a d-wave order parameter in the system studied, namely 
YBa2Cu307_,5 with small 5. Spin-fluctuation models of high-T c lead to d-wave supercon- 
ductivity Jl||. Strong correlation models have a similar instability to d-wave pairing H0. 
Thus the experiments on YB^CusCv-^ may be taken as a confirmation of the magnetic 
mechanism in general. To distinguish between the two models is not an easy problem. 
However, variational Monte Carlo (VMC) calculations on a strong correlation model, the 
extended t-J model, have suggested that, while a d-wave is the dominant instability, s-wave 
and d-wave mixing may be the favored ground state for higher doping levels || and certain 
band structure parameters ||. From this point of view, not all high-T c systems are alike, 
and it becomes crucial to perform experiments on different systems and, most importantly, 
at different doping levels. 

Apart from the tunneling technique, the most direct measurement (although not phase- 
sensitive) is angle-resolved photoemission (ARPES). A recent experiment by Ma et al. JTO 



used ARPES to investigate the gap in Bi 2 Sr 2 CaCu208+z single crystals. In particular this 
group has measured the gap at two points on the Fermi surface: A(T — M) or A(k x = 0) in 
our coordinate system (along the Cu-0 bond) and A(T — X) or A(k x = k y ). This was done 
as a function of temperature. The basic features of the experimental results are : (i) the 
gap A(r — M) opens up quickly below T c and A(r — M) > A(T — X) for all temperatures; 
(ii) A(r — X) is small or even zero for T > 0.8 T c and then increases gradually as T is 
further reduced, leading to a strongly temperature-dependent gap anisotropy; (iii) there are 
indications of upward curvature in A{Y—X); (iv) the temperature dependence of A is weak in 
both directions for a wide range of temperatures (from T=0 K to T=0.7 T c approximately). 



These data are different from results obtained by other workers on B^S^CaC^Og+a; JTT 



in that the gap along r — X is not always small. This is due to annealing in an oxygen 



atmosphere, leading to higher oxygen content and the presumably to a higher hole doping 
level in the Cu0 2 planes. [K| 

We do not possess a detailed microscopic theory of s-d mixing and the VMC calculations 
are restricted to T=0. Thus we work within the framework of Ginzburg-Landau theory. The 
point group symmetry of the CuO"2 planes in the temperature range of interest is orthorhom- 
bic C2V1 due to buckling of the planes fl3] . In this symmetry the interesting thing is that an 



order parameter with d-wave symmetry transforms according to the identity representation. 
[These symmetry issues are discussed in detail in Ref. [|J]. So from a strict symmetry point 
of view a d-wave and an s-wave are indistinguishable. The concept of a two-component 
order parameter, where one component has an s-wave symmetry and the other a d-wave 
symmetry is thus well-defined only in the tetragonal case. In the present orthorhombic case 
we can regard the s-d mixing as a strong temperature dependent anisotropy which neverthe- 
less never actually breaks the crystalline symmetry. If we write the Ginzburg-Landau free 
energy density as follows: 

/ = a° s (T - T S )\V S \ 2 + a° d (T - T d )\^ d \ 2 + (3 s \^ s f + (3 d \V d \ 4 + (3\^ s \ 2 \^ d \ 2 

+(3 sd (y s 2 V d * 2 + c.c.) + 7 + c.c.) , (1) 

the difference from the case of tetragonal symmetry is the addition of the final bilinear term. 
In fact the coupling constant 7 is proportional to the orthorhombic distortion b — a where b 
and a are the lattice constants in the Cu-0 plane. Note that in our case, the gradient terms 
are not necessary due to the absence of external fields and the assumed homogeneity of the 
material (translational invariance). We minimize the free energy with respect to s \ and 
\^d\ arid the phase difference of the two components 0. The phase difference is crucially 
important. If it is zero then the resulting state is s+d: the gap function is real (except for 
the overall phase related to gauge freedom). If <ft is nonzero the state is more similar to s+id: 
it breaks time reversal symmetry and has no nodes. 

The minimization leads to coupled equations which we have solved numerically. The 
result of the calculations is that, according to the sign of the parameter (3 s d in the above 



functional, there are two possibilities for the value of 0. It is either zero or there is a 
transition from a continuous value (temperature dependent) to zero. We attempted to fit 
the data with a temperature dependent phase shift without success. The data therefore 
point to the s+d state. So we adopted the case = for the phase difference between the 
two components of the order parameter. We chose (3 sd < (which then essentially adds to 
/3) and 7 < 0. The equations we have to solve in order to get the order parameters become: 

a s>d \^ s>d \ + 2/3 M |# S)d | 3 + (J3 + (3 sd ) \V d , s \ 2 \y s4 \ + ||* d , s | = 0. (2) 

The gap function we consider is of the form: 

A(k) = |* d | f d (k) + / s (k), (3) 

where k lies on the Fermi surface. The assumed forms of the k-dependence are : 

f d (k x = k y ) = 0J d (k x = 0,k y ) = l (4) 
fs(k x = ky) = f s (k x = 0, ky) = 1 (5) 

Accordingly, the gap along the T—X direction is the solution \ ^> s \ whereas along the T—M 
direction it is |\I/ S | + l^l- In order to compare the data over the whole temperature range 
to experiment we used the approximation a S)d / '/3 S)d oc [tanh(z/(^ — 1) 1//2 )] 2 . For v = 1.74 
this form gives a good fit to the weak-coupling BCS temperature dependence of the gap. 
Increasing v is a way of mocking up strong-coupling effects. In Fig. [I] we show the results 
of the fitting to experimental data. Several features of the fit are noteworthy. A strong- 
coupling value of the parameter v must be used in order to represent feature (i) mentioned 
above. T s must be lower than T^ by about 15 K so that the s-component is suppressed near 
T c , feature (ii) above. The upward curvature in A(T — X) [feature (iii)] is only produced 
by having a nonzero value of 7, the bilinear coupling of the two components. Finally, the 
tangent function in the coefficients of the quadratic terms in the free energy produces a 
temperature independent gaps at low T, as in ordinary BCS theory, [feature (iv)]. There is 
only one second order transition. 
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Near T c the gap is pure d-wave, with gap zeros in the diagonals and maxima along 
the axes. As the temperature is reduced, the s-wave component grows, the nodes move 
towards the k y axis and the gap along this axis become smaller while the gap along the 
k x axis becomes larger. Thus it is a temperature-dependent orthorhombic anisotropy. We 
take, for illustration purposes only, a circular Fermi surface with f d = cos(29k) , f s = 1. 
Then the zeros of the gap move from a value of 6 k = 7r/4 close to T c to a value which is 
temperature dependent ( 9k = — 1/2 arccos(— ^ s /^d) )• This is shown in Fig. |2|. In this 
extreme example the zeros of the gap move all the way to the k y axis and then disappear. 
We stress that this disappearance is a highly model-dependent phenomenon, depending, for 
example on the detailed momentum dependence of the s-wave component. On the other 
hand, the qualitative result that the zeros of the gap will move as the temperature is varied 
is a firm prediction of the two-component theory as long as the components have different 
temperature dependences, as is indicated by the experiments. Note also that experiments 
such as the reported in Ref.l which depend only on the relative phase change under a n/2 
rotation in k-space and will not be affected by a small movement of the nodes. 

In Fig. [3] , the anisotropy ratio A(T — X)/A(T — M) is plotted as a function of tem- 
perature. The dramatic temperature variation of this quantity is well reproduced by the 
strong-coupling parameters. Also in Fig. |], the form of the heat capacity (in arbitrary units) 
is plotted. Since we are interested for values close to T c , we used the linear form of the 
parameters a s and , in such a way that the first term of the Taylor expansion (with T c 
instead of T in the denominator) of the hyperbolic tangent agrees. The hyperbolic tangent 
form was chosen above in order to compare the data with the theory at much lower tem- 
peratures. A somewhat anomalous behaviour of C V (T) is predicted, with a local maximum 
around 0.7 T c . This comes from the sharp growth of the second component in the order 
parameter in this temperature range. 

We do not claim that the fit in Fig. [I] is a quantitative confirmation of the two-component 
theory. The curves are smooth and the numbers of parameters large. It is true, however, 
that the temperature-dependent anisotropy of the superconducting gap in BSCCO is what 
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one expects if a two-component order parameter exists. A pure d-wave order parameter 
alone cannot account for the data. Our main point in this paper is to note that the two- 
component picture has a number of surprising consequences which may be checked in highly 
hole-doped BSCCO. These are : the gap energy should have orthorhombic anisotropy at low 
temperatures; there should be upward curvature in A(T — X)(T) already suggested by the 
existing data; the position of the gap nodes should be temperature-dependent; there should 
be a peak in the electronic specific heat close to 0.7 T c . 

We would like to thank Marshall Onellion, Ron Kelley, Jian Ma, Christoph Quitmann for 
valuable discussions. This work was supported by the National Science Foundation under 
grants No. 9214739 and No. 921047. 
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FIGURES 

FIG. 1. Shows the gap along the two directions of the Fermi surface in comparison to the 
experimental data. The solid line is the solution for v = 3.0 (the values of the parameters 
are : a s /p s = -250[tanh(3.0(0.8T c /T - l) 1 / 2 )] 2 , 7 /ft = 7 /Ai = "100, P/P s = P/P d = 2, 
a d /p d = -220[tanh(3.0(T c /T - 1) 1/2 )] 2 , the dashed line is the solution for v = 1.74 (the val- 
ues of the parameters are: a s /(3 s = -340[tanh(1.74(0.8T c /T - 1) 1/2 )] 2 , j/P s = j/(3 d = -100, 
P/Ps = P/Pd = 2, a d /p d = -280[tanh(1.74(T c /T - I) 1 / 2 )} 2 ). 

FIG. 2. Shows the zeros of the gap as a function of temperature, for the hypothetical circular 
Fermi surface. At low temperatures there are no zeros in the gap according to the model and ^ s 
is greater than ty d . 

FIG. 3. Shows the anisotropy ratio defined as A(r - X)/A(T - M) for the value of v = 3.00. 
FIG. 4. The specific heat as predicted by the model. Near 0.7 T c there is a local maximum. 
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